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Exploring positron characteristics utilizing two new positron-electron correlation schemes based on 

multiple electronic-structure calculation methods 
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We make a gradient correction to a new local density approximation form of positron-electron correlation. 

Then the positron lifetimes and affinities are probed by using these two approximation forms based on three 
electronic-structure calculation methods including the full-potential linearized augmented plane wave (FLAPW) 
plus local orbitals approach, the atomic superposition (ATSUP) approach and the projector augmented wave 
(PAW) approach. The differences between calculated lifetimes using the FLAPW and ATSUP methods are 
clearly interpreted in the view of positron and electron transfers. We further find that a well implemented PAW 
method can give near-perfect agreement on both the positron lifetimes and affinities with the FLAPW method, 
and the competitiveness of the ATSUP method against the FLAPW/PAW method is reduced within the best 
calculations. By comparing with experimental data, the new introduced gradient corrected correlation form is 
proved competitive for positron lifetime and affinity calculations. 

PACS numbers: 78.70.Bj, 71.60.+Z, 71.15.Mb 


I. INTRODUCTION 

In recent decades, the Positron Annihilation Spectroscopy 
(PAS) has become a valuable method to study the microscopic 
structure of solids and gives detailed information on the 
electron density and/or momentum distribution |3l in the re¬ 
gions scanned by positrons. An accompanying theory is re¬ 
quired for a thorough understanding of experimental results. 
A full two-component self-consistent scheme W has been 
developed for calculating positron states in solids based on 
the density functional theory (DFT) ||3l. Especially in bulk 
material where the positron is delocalized and does not af¬ 
fect the electron states, the full two-component scheme can be 
reduced without losing accuracy to the conventional scheme 
W in which the electronic-structure is determined by com¬ 
mon one-component formalism. However, there are various 
kinds of approximations can be adjusted within this calcu¬ 
lations. To improve the analyses of experimental data, one 
should find out which approximations are more credible to 
produce the positron state llsl- fl^ . In this short paper, we fo¬ 
cus on probing the positron lifetimes and affinities by using 
two new positron-electron correlation schemes based on three 
electronic-structure calculation methods. 

Recently, N. D. Drummond et al. made two cal¬ 

culations for a positron immersed in a homogeneous elec¬ 
tron gas, by using the Quantum Monte Carlo (QMC) method 
and a modified one-component DFT method, and then two 
forms of local density approximations (LDA) on the positron- 
electron correlation are derived. Kuriplach and Barbiellini 
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JH lil] proposed a fitted LDA form and a generalized gradi¬ 
ent approximation (GGA) form based on previous QMC cal¬ 
culation, and then applied these two forms to multiple cal¬ 
culations for positron characteristics in solid. However, the 
LDA form based on the modified one-component DFT calcu¬ 
lation has not been studied. In this work, we make a gradi¬ 
ent correction to the IDFTLDA form and validate these two 
new positron-electron correlation schemes by applying them 
to multiple positron lifetimes and affinities calculations. 

Besides, we probe in detail the effect of different electronic- 
structure calculation methods on positron characteristics in 
solid. These methods include the full-potential linearized aug¬ 
mented plane-wave (FLAPW) plus local orbitals method IlLJIl . 
the projector augmented wave (PAW) method ifldll . and the 
atomic superposition (ATSUP) method iflsll . Among these 
methods, the FLAPW method is regarded as the most ac¬ 
curate method to calculate electronic-structure, the ATSUP 
method performs with the best computational efficiency, the 
PAW method has greater computational efficiency and close 
accuracy as the FLAPW method but has not been com¬ 
pletely tested on positron lifetimes and affinities calcula¬ 
tions except some individual calculations iflhUl^ . More¬ 
over, our previous work 1^ showed that the calculated 
lifetimes utilizing the PAW method disagree with that uti¬ 
lizing the FLAPW method. However, within those PAW 
calculations, the ionic potential was not well constructed. 
In this paper, we investigated the influences of the ionic 
pseudo-potential/full-potential and different electron-electron 
exchange-correlations approaches within the PAW calcula¬ 
tions. Especially, the difference between calculated life¬ 
times by using the self-consistent (ELAPW) and non-self- 
consistent (ATSUP) methods is clearly investigated in the 
view of positron and electron transfers. 

This paper is organized as follows: In Sec. 2, we give a 
brief and overall description of the models considered here 
as well as the computational details and the analysis meth¬ 
ods we used. In Sec. 3, we introduce the experimental data 
on positron lifetime used in this work. In Sec. 4, we firstly 
apply all approximation methods for electronic-structure and 
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positron-state calculations to the cases of Si and Al, and 
give detailed analyses on the effects of these different ap¬ 
proaches, and then assess the two new correlation schemes 
by using the positron lifetime/affinity data in comparison with 
other schemes base on different electronic-structure calcula¬ 
tion methods. 


II. THEORY AND METHODOLOGY 


A. Theory 


In this section, we briefly introduce the calculation scheme 
for the positron state and various appproximations investi¬ 
gated in this work. Firstly, we do the electronic-structure cal¬ 
culation without considering the perturbation by positron to 
obtain the ground-state electronic density ne^{f) and Coulomb 
potential Vcouii^ sensed by positron. Then, the positron den¬ 
sity is determined by solving the Kohn-Sham Eq.: 

[-]ff+Vc„ul{f) + Vcorr(.^]ilf^ = ( 1 ) 

where Vcorri^ is the correlation potential between electron 
and positron. Finally, the positron lifetime can be obtained 
by the inverse of the annihilation rate, which is proportional 
to the product of positron density and electron density accom¬ 
panied by the so-called enhancement factor arising from the 
correlation energy between a positron and electrons 1^ . The 
equations are written as follows: 


j, A^nr^cJdrne-{i^ne+(f)7{ne-), 


( 2 ) 


where is the classical electron radius, c is the speed of light, 
and is the enhancement factor of the electron density 

at the position f. The positron affinity can be calculated by 
adding electron and positron chemical potentials together: 


A"* = H-P+. 


(3) 


The positron chemical potential is determined by the 
positron ground-state energy. The electron chemical poten¬ 
tial fT is derived from the Fermi energy (top energy of the 
valence band) in the case of a metal (a semiconductor). This 
scheme is still accurate for a perfect lattice, as in this case the 
positron density is delocalized and vanishingly small at every 
point thus does not affect the bulk electronic-structure 1^ 1^ . 


TABLE I: Parameterized LDA/GGA correlation schemes. 


y 

ai 

(33 

03/2 

07/3 

08/3 

a 

IDFTLDA 

4.1698 

0.1737 

-1.567 

-3.579 

0.8364 

0 

IDFTGGA 

4.1698 

0.1737 

-1.567 

-3.579 

0.8364 

0.143 

fQMCLDA 

-0.22 

1/6 

0 

0 

0 

0 

fQMCGGA 

-0.22 

1/6 

0 

0 

0 

0.05 

PHCLDA 

-0.137 

1/6 

0 

0 

0 

0 

PHCGGA 

-0.137 

1/6 

0 

0 

0 

0.10 


In practice of this work, each enhancement factor is applied 
identically to all electrons as suggested by K. O. Jensen 0. 
These enhancement factors can be divided into two categories: 
the local density approximation (FDA) and the generalized 
gradient approximation (GGA), and parameterized by the fol¬ 
lowing equation, 

y= \ + {\.22>rs +a2r] + a^r] + aii2r]^^ 

+a2i2,rj^ + (4) 

here, is defined by - (3/47rne,_)'''^, e is defined by 
e = |V ln(ne,_)|^/q'^p (q'^p is the local Thomas-Fermi screening 
length), a 2 , ag, ag/a. asji, ^7/3. ag/a and a are fitted parame¬ 
ters. We investigated five forms of the enhancement factor and 
correlation potential marked by IDFTLDA ifl^ . fQMCLDA 
HH, fQMCGGA Hi], PHCLDA 0 and PHCGGA 0, 
plus a new GGA form IDFTGGA introduced in this work 
based on the IDFTLDA scheme. The fitted parameters of 
these enhancement factors are listed in Table U . The FDA 
forms of Vcorr correspondin g to IDFTLDA, f Q MCLDA, PH¬ 
CLDA are given in Refs. llT3l . H and ll2-5n . respectively. 
Within the GGA, the correspon ding correlation potential takes 
the form The electronic den¬ 

sity and Coulomb potential were calculated by using various 
methods including: a) the all-electron full potential linearized 
alimented plane wave plus local orbitals (FLAPW) method 
ifo as implemented in Ref.H being regarded as the most 
accurate method to calculate electronic-structure, b) the pro¬ 
jector augmented wave (PAW) method ifldl] with reconstruc¬ 
tion of all-electron and full-potential performing with greater 
computational efficiency and close accuracy as the FLAPW 
method, c) the non-self-consistent atomic superposition (AT- 
SUP) method ifl^ performing with the best computational ef¬ 
ficiency. 


B. Computational details 

During the calculations for electronic-structure, three meth¬ 
ods mentioned above are implemented in this work. For 
FLAPW calculations, the WIEN2k code 10 was used, the 
PBE-GGA approach 10 was adopted for electron-electron 
exchange-correlations, the total number of k-points in the 
whole Brillouin zone (BZ) was set to 3375, and the self- 
consistency was achieved up to both levels of 0.0001 Ry for 
total energy and 0.001 e for charge distance. For PAW cal¬ 
culations, the PWSCF code within the Quantum ESPRESSO 
pack age ll^ was used, the PBEsol-GGA 1I3II1 and PZ-LDA 
113^ approaches were also implemented for electron-electron 
exchange-correlations besides the PBE-GGA approach, the 
PAW pseudo-potential files named PSLibrary 0.3.1 and gen¬ 
erated by A. D. Corso (SISSA, Italy) were employed 
the k-points grid was automatically generated with the param¬ 
eter being set at least (333) in Monkhorst-Pack scheme, the 
kinetic energy cut-off of more than 100 Ry (400 Ry) for the 
wave-functions (charge density) and the default convergence 
threshold of 10 ® were adopted for self-consistency. For AT- 
SUP calculations, the electron density and Coulomb potential 
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for each material were simply approximated by the superposi¬ 
tion of the electron density and Coulomb potential of neutral 
free atoms Q, while the total number of the node points 
was set to the same as in PAW calculations. Besides, the 
2x2x2 supercells were used to calculate the electron struc¬ 
tures of monovacancy in A1 and Si. To obtain the positron- 
state, the three-dimensional Kohn-Sham equation Eq. ([TJ was 
solved by the finite-difference method while the unit cell of 
each material was divided into about 10 mesh spaces per bohr 
in each dimension. All important variable parameters were 
checked carefully to achieve that the computational precision 
of lifetime and affinities are the order of 0.1 ps and 0.0 leV, 
respectively. 

C. Model comparison 

An appropriate criterion must be chosen to make a com¬ 
parison between different models. The root mean squared de¬ 
viation (RMSD) is the most popular one and defined as the 
square root of the mean of the squared deviation between ex¬ 
perimental and theoretical results: RMSD = “ 

^theo^ 2 ^^j 1/2, here N denotes the number of experimental val¬ 
ues. In addition, since the theoretical values can be treated to 
be noise-free, the simple mean-absolute-deviation (MAD) de¬ 
fined by MAD = is much more 

meaningful to quantify the overall differences between calcu¬ 
lated results by using various models. It is obvious that the 
experimental data favor models producing lower values of the 
RMSD. 


III. EXPERIMENTAL DATA 

Up to five recent observed values from different literatures 
and groups for 21 materials were gathered to compose a re¬ 
liable experimental data set. All the experimental values for 
each material investigated in this work are collected basically 
by using the standard suggested in Ref. ifs^ and listed in 
Table HJ with their standard deviation. Furthermore, the mate¬ 
rials having less than five experimental measurements and/or 
the older experimental data were avoided being adopted. It 
is reasonable to suppose that these materials having insuffi¬ 
cient and/or unreliable experimental data would disrupt the 
comparison between models. Especially, the measurements 
for alkali-metals reported before 1975 are not suggested to 
be treated seriously |@]. The deviations of experimental re¬ 
sults between different groups are usually much larger than 
the statistical errors, even when just the recent and reliable 
measurements are considered. That is, the systematic error is 
the dominant factor, so that the sole statistical error is far from 
enough and not used in this work. However, the systematic er¬ 
ror is difficult to derive from single experimental result. So in 
this paper, the average experimental values of each material 
were used to assess the positron-electron correlation models, 
and the systematic errors are expected to be cancelled as in 
Ref. Js^l- Because the observed values for defect state are 
insufficient and/or largely scattered, it is hard to make a clear 


TABLE II: The experimental values of lifetime Texp, the related mean 
value T*xp and the corresponding standard deviation CTexp for each 
material involved in this work. 


Material 

^exp 

^exp 

^exp 

Si 

216.7[^ 218[3A| 218[^ 222[^ 216[3A| 

218.1 

2.323 

Ge 

220.5\3£ 230[^ 230[3£| 228[3£1 228[3£| 

227.3 

3.931 

Mg 

225[^ 225[3^ 220[^ 238[^ 235[34] 

228.6 

7.569 

A1 

160.7[3£1 166[^ 163[3L| 165[3£1 165[^ 

163.9 

2.114 

Ti 

147[^ 154[^ 145[3L| 152[3£] 143[^ 

148.2 

4.658 

Fe 

108[^ 106[3^ 114[^ 1I0[^ ni[3^ 

109.8 

3.033 

Ni 

109.8[3£1 107[3^ I05|3£] 109[3£] 110[3^ 

108.2 

2.127 

Zn 

148[3^ 153|^ 145[3A| 154[3£] 1521^ 

150.4 

3.781 

Cu 

110.7[^ 122\3^ 112[3L| 110[^ 120[^ 

114.9 

2.514 

Nb 

119[^ 120[3A| I22[3L| 122[^ 125[^ 

121.6 

2.302 

Mo 

109.5[^ 103[^ 118[3A| 114[^ 104[34] 

109.7 

6.418 

Ta 

1161^ 122[3^ 120[3L| 1251^ 117[34] 

120.0 

3.674 

Ag 

120|3£1 1301^ 131[3L| 133|^ 131|^ 

129.0 

5.147 

Au 

1I7[3£] 1131MI I13[3L] 1I7[3£] 123[3£| 

116.6 

4.098 

Cd 

175[^ 184[3A| I67[^ 172[^ 186[3^ 

176.8 

8.043 

In 

194.7[^ 200[3£| 192[3L| 193[^ 189[^ 

193.7 

4.066 

Pb 

194[^ 200[3^ 204[^ 200[^ 209[^ 

201.4 

5.550 

GaAs 

231.6[3T] 23 1[^ 230[3£| 232[^ 220[41J 

228.9 

5.043 

InP 

2411^ 240[^ 247[^ 2421^ 244[^ 

242.8 

2.775 

ZnO 

153|^ 159[^ 158[^ 161[^ 1711^ 

160.4 

6.618 

CdTe 

284[^ 2851^ 285[^ 289[«] 29I[56] 

286.8 

3.033 


discussion on the defect state by using these positron-electron 
correlation models in this short paper. Thus, except the de¬ 
tailed analyses in the cases of Si and A1 based on three usu¬ 
ally applied approaches for electronic-structure calculations, 
we mainly focus on testing the correlation models by using 
bulk materials’ lifetime data and positron-affinity data. The 
experimental data of positron affinity are listed in Table IV] 


IV. RESULTS AND DISCUSSION 
A. Detailed analyses in cases of Si and A1 

Representatively, the panels (a) and (c) in Fig. [T] (Fig. |3) 
show respectively the self-consistent all-electron and positron 
densities on plane (110) for A1 (Si) based on the FFAPW 
method together with the fQMCGGA form of the enhance¬ 
ment factor and correlation potential. It is reasonable to ob¬ 
tain that the panel (a) in Fig. |2] shows clear bonding states 
of Si while the panel (a) in Fig. [T] shows the presence of 
the nearly free conduction electrons in interstitial regions. To 
make a comparison between the FFAPW and ATSUP method 
for electronic-structure calculations, we also plot the ratio of 
their respective all-electron and positron densities in panel (b) 
and (d) in Fig. [T|(Fig. |2]l for A1 (Si). These four ratio panels 
actually reflect the electron and positron transfers from densi¬ 
ties based on the non-self-consistent free atomic calculations 
to that based on the exact self-consistent calculations. It con¬ 
firms the fact that the positron density follows the changes of 
the electron density which yield not a big difference in anni- 
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FIG. 1; Left panels: the self-consistent all-electron (a) and positron 
densities (c) (in unit of a.u.) on plane (110) for A1 based on the 
FLAPW method and the fQMCGGA approximation. Right panels: 
the ratios of all-electron (b) or positron densities (d) calculated by 
using the FLAPW method to that by using the ATSUP method. 



[100] (Angstrom) 



[100] (Angstrom) 


0) 

Q. 



FIG. 2: AsFig.[I]butforSi. 


hilation rate between these two calculations 

Now, taking more subtle analyses, the change of lifetime 
within the FLAPW calculation from that within the ATSUP 
calculation for A1 is attributed to the competition between the 
following two factors: a) the lifetime is decreased by the trans¬ 
lations of electrons (illustrated in Fig. [1] (b) as T^) from 
near-nucleus regions with tiny positron densities to intersti¬ 
tial regions with large positron densities, b) the lifetime is in¬ 
creased by the translation of positron (illustrated in Fig. [T|(d) 
as from core regions with large electron densities to in¬ 
terstitial regions with small electron densities. However, in 
the case of Si with bonding states, the change of lifetime de¬ 
pends conversely on the translations of electrons and positron: 
a) the lifetime is increased by the translations of electrons (il¬ 
lustrated in Fig. |2](b) as from interstitial regions with the 
largest positron densities to bonding regions with tiny positron 


FIG. 3: The total Coulomb potential (in unit of Ry) sensed by 
the positron based on the ionic pseudo-potentials (Vpp) and recon¬ 
structed ionic full-potential (Vpp) and the corresponding calculated 
positron densities (in unit of a.u.) along the [100] direction be¬ 
tween two adjacent atoms for A1 (a) and Si (b), respectively. To 
make a further comparison, the full-potentials calculated by using 
the FLAPW method (VppApw) are also plotted. 


densities, b) the lifetime is decreased by the translation of 
positron (illustrated in Fig. |2] (d) as T^p from interstitial 
regions with tiny electron densities to bonding regions with 
large electron densities. Taking note of the magnitude of scale 
rulers, these two figures state clearly that the translations of 
electrons (T*^') are dominant factors for both A1 and Si. Con¬ 
sequently, the lifetimes within the FLAPW calculations be¬ 
come smaller (larger) for A1 (Si). These variances are proved 
by calculated values of lifetimes listed in Table [HI] In addi¬ 
tion, the lifetimes of Si calculated by using three GGA forms 
of the enhancement factor show greater differences since the 
large electron-density gradient terms in bonding regions giv¬ 
ing decreases of the enhancement factor can further weaken 
the effect of the translation T^p 

We calculated the bulk lifetimes for A1 and Si based on the 
PAW method. Within the Table HIH the label ”PAW” without 
a suffix indicates that the electron-structure is calculated by 
using the PBE-GGA electron-electron exchange-correlations 
approach ll2^ and positron-state is calculated by using re¬ 
constructed ionic full-potential (FP), the suffix ”-PZ” indi¬ 
cates that the PBE-GGA approach is replaced by the PZ-LDA 
approach during electron-structure calculations, and the 
sufhx ”-PP” indicates that ionic full-potential (EP) is replaced 
by the ionic pseudo-potential (PP) during positron-state cal¬ 
culations. The ionic potential together with the Hartree po¬ 
tential from the valence electrons compose the total Coulomb 
potential in Eq. ([T]i. It can be easily found that the better im¬ 
plemented PAW method by using reconstructed full-potential 



















5 


can give a startling agreement with the FLAPW method on the 
positron-lifetime calculations for A1 and Si. By comparing the 
results of PAW and PAW-PP approach, the PAW-PP approach 
leads to smaller lifetimes with the differences up to 3.8 ps and 
4.3 ps for A1 and Si respectively. These decreases are caused 
by the fact that the softer potential within the PAW-PP ap¬ 
proach more powerfully attracts positron into the near-nucleus 
regions with much larger electron densities. This statement is 
illustrated by the Fig. [3 showing the total Coulomb potential 
sensed by the positron based on the ionic pseudo-potential 
(Vpp) and reconstructed ionic full-potential (Vpp) and the cor¬ 
responding calculated positron densities p®^ along the [100] 
direction between two adjacent atoms for A1 (a) and Si (b), re¬ 
spectively. To make a further comparison, the full-potentials 
calculated by using the FLAPW method (Vflapw) are also 
plotted and found nearly the same as the reconstructed PAW 
full-potentials. This figure indicates that a change in the ionic 
potential approachs (FP or PP) can lead to a change of more 
than one order of magnitude in the positron densities near 
the nuclei. It should be noted that, in cases of PAW calcu¬ 
lations with underestimated core/semicore electron densites 
in the near-nucleus regions ifs^ . the effect of overestimated 
positron densities based on the pseudo-potentials can be can¬ 
celled, and then excellent quality on the calculated positron 
lifetimes is able to be achieved. It is clear that the differences 
between the results of PAW-PZ and PAW are of the order of 
0.1 ps, and therefore the effect of different electron-electron 
exchange-correlations schemes is small. More than this, we 
also calculated the lifetimes by using the PBEsol-GGA ap¬ 
proach M which is revised for solids and their surfaces, and 
the similar differences of the order of 0.1 ps are also obtained 
compared with the PBE-GGA approach. 


TABLE III: Calculated results of positron lifetimes (in unit of ps) for 
Al, Si, and ideal monovacancy in A1 and Si based on various methods 
for electronic-structure and positron-state calculations. 




IDFT 

IDFT 

fQMC 

fQMC 

PHC 

PHC 



GGA 

LDA 

GGA 

LDA 

GGA 

LDA 


ATSUP 

160.778 

152.470 

173.347 

169.357 

163.036 

156.438 

Al 

FLAPW 

156.615 

149.852 

169.972 

166.530 

159.397 

153.878 

PAW 

156.649 

149.898 

170.016 

166.584 

159.432 

153.925 


PAW-PP 

154.113 

146.814 

166.507 

162.798 

156.574 

150.587 


PAW-PZ 

157.208 

150.204 

170.421 

166.906 

159.898 

154.220 


ATSUP 

201.770 

186.634 

213.260 

207.345 

201.363 

190.484 

Si 

FLAPW 

211.843 

188.285 

217.520 

208.477 

208.639 

191.790 

PAW 

211.779 

188.245 

217.466 

208.431 

208.586 

191.752 


PAW-PP 

208.407 

184.675 

213.320 

204.125 

205.060 

187.976 


PAW-PZ 

211.248 

188.388 

217.399 

208.625 

208.247 

191.905 

Vai 

ATSUP 

229.441 

216.639 

246.294 

240.941 

229.686 

220.274 

PAW 

212.176 

201.245 

229.481 

224.429 

214.050 

205.570 

Vsi 

ATSUP 

227.458 

208.972 

239.524 

232.309 

225.922 

212.690 

PAW 

236.052 

208.712 

241.816 

231.443 

231.504 

212.145 


In addition, as shown in Table [111] the positron lifetimes 
for monovacancy in Al and Si are also calculated based on 
the ATSUP and PAW methods for electronic-structure calcu¬ 
lations and six correlation schemes for positron-state calcula¬ 


tions. The ideal monovacancy structure is used in these cal¬ 
culations, which means that the positron is trapped into a sin¬ 
gle vacancy without considering the ionic relaxation from the 
ideal lattice positions. Larger differences between the results 
of ATSUP and PAW are found in monovacancy-state calcula¬ 
tions compared with that in bulk-state calculations. Besides, 
the IDETGGA/IDETLDA correlation schemes produce sim¬ 
ilar lifetime values compared with the PHCGGA/PHCLDA 
correlation schemes and produce much smaller lifetime val¬ 
ues compared with the fQMCGGA/fQMCLDA correlation 
schemes in both monovacancy-state and bulk-state calcula¬ 
tions. 


B. Positron lifetime calculations 

In this subsection we firstly give visualized comparisons be¬ 
tween experimental values and calculated results based on dif¬ 
ferent methods for electronic-structure and positron-state cal¬ 
culations. Within the PAW, the positron lifetimes are all cal¬ 
culated by using the reconstructed full-potential and certainly 
all-electron densities from now on. 



FIG. 4: The deviations of the theoretical results based on various 
methods from the experimental values alongwith the standard devia¬ 
tion of experimental values for each material. 


TABLE IV: The MADs between the calculated results by using the 
ATSUP/PAW method and that by using the FLAPW method. And 
the RMSDs between the theoretical results and the experimental data 

Axp* 



MAD [ps] 

ATSUP PAW 

FLAPW 

RMSD [ps] 

PAW 

ATSUP 

fQMCGGA 

2.503 

0.303 

4.503 

4.591 

6.309 

IDFTGGA 

5.068 

0.316 

4.809 

4.821 

5.611 

PHCGGA 

3.667 

0.287 

6.148 

6.013 

7.672 

fQMCLDA 

2.184 

0.290 

11.36 

11.19 

10.35 

IDFTLDA 

1.966 

0.253 

25.19 

24.99 

23.88 

PHCLDA 

1.936 

0.260 

22.83 

22.63 

21.54 
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The deviations of the theoretical results from the experi¬ 
mental data alongwith the standard deviations of observed val¬ 
ues for all materials are plotted in Fig. |4] The scattering re¬ 
gions of calculated results by different forms of the enhance¬ 
ment factor are found much larger in the atom systems with 
bonding states compared with that in pure metal systems. Be¬ 
sides, the deviations of the results by using the ATSUP method 
from those by using the FLAPW method are mostly larger in 
GGA approximations compared with those in LDA approxi¬ 
mations. Numerically, the MADs for different forms of the 
enhancement factor between the calculated lifetimes by using 
the ATSUP method and those by using the FLAPW method 
are shown in Table |IV] These MADs range from 1.936 ps 
(PHCLDA) to 5.068 ps (IDFTGGA). Moreover, the well im¬ 
plemented PAW method is found able to give nearly the same 
results as the FLAPW method. Numerically, the MADs be¬ 
tween the calculated lifetimes by the PAW method and those 
by the FLAPW method for different forms of the enhancement 
factor are also shown in Table |IV] These MADs range from 
0.253 ps (IDFTLDA) to 0.316 ps (IDFTGGA). This near¬ 
perfect agreement between the PAW method and the FLAPW 
method proves our calculations are quite credible. 

Table IIVI also presents the RMSDs between the theoretical 
results and the experimental data T*^p by using six positron- 
electron correlation schemes. Two interesting phenomena 
can be found in this table. Firstly, the RMSDs produced by 
the IDFTLDA scheme are always worse among the RMSDs 
based on three electron structure approaches, but are similar 
to those produced by the PHCLDA scheme. Thus, the gra¬ 
dient correction (IDFTGGA) to this LDA form (IDFTLDA) 
is needed. It is clear that the corrected IDFTGGA scheme 
largely improves the calculations, and performs better than 
the PHCGGA scheme, but is still worse than the fQMCGGA 
scheme. The fQMCGGA scheme together with the FLAPW 
method produced the best RMSD. This fact indicates that the 
quantum Monte Carlo calculation implemented in Ref. iflTl] 
is more credible than the modified one-component DFT cal¬ 
culation lfT3] on the positron-electron correlation. Secondly, 
compared to the RMSD produced by using the FLAPW/PAW 
method, the RMSD produced by using the simple ATSUP 
method is a little smaller based on the LDA correlation 
schemes, but is distinctly larger based on the GGA (espe¬ 
cially fQMCGGA) correlation schemes. This phenomenon 
implies that the benefit of the exact eletronic-structure cal¬ 
culation approach (PAW/FLAPW) is swamped by the in¬ 
accurate approximation of the enhancement factor. Mean¬ 
while, the competitiveness of the ATSUP approach against the 
FLAPW/PAW method is reduced based on the most accurate 
positron-electron correlation schemes. 


C. Positron affinity calculations 

The positron affinity is a important bulk property which 
describes the positron energy level in the solid, and allows us 
to probe the positron behavior in an inhomogeneous material. 
For example, the difference of the lowest positron energies be¬ 
tween two elemental metals in contact is given by the positron 


affinity difference, and determines how the positron samples 
near the interface region. Besides, if the electron work func¬ 
tion </>“ is known, the positron work function can be de¬ 
rived by the equation: . The crystal (e.g., W 

metal) having a stronger negative positron work function can 
emit slow-positron to the vacuum from the surface and there¬ 
fore be utilized as a more efficient positron moderator for the 
slow-positron beam. 

The theoretical and experimental positron affinities for 
eight common materials by using the new IDFTLDA and 
IDFTGGA correlation schemes are listed in Table FV] To make 
a comparison, the results corresponding to the PHCGGA and 
fQCMGGA schemes are also listed. During the electron struc¬ 
ture calculation, the ATSUP method was not implemented 
because the ATSUP method is inappropriate for positron en¬ 
ergetics calculations and gives much negative positron work 
functions ifl^ . Within the PAW calculations, both the PBE- 
GGA and PZ-LDA approaches are used for electron-electron 
exchange-correlations. The RMSDs between the theoretical 
and experimental positron affinities are also presented in Ta¬ 
ble |V] 

As in previous lifetime calculations, the calculated positron 
affinities by using the FLAPW method are also near the same 
as that by using the PAW method. Besides, our calculated 
positron affinities by using the fQMCGGA & PZ-LDA ap- 
OToaches are in excellent agreement with that reported in Ref. 
IK] with a MAD being 0.06 eV. Moreover, the differences 
between the RMSDs produced by using the PBE-GGA and 
PZ-LDA approaches, are not negligible, and the PBE-GGA 
approach performs mostly better than the PZ-LDA approach 
except the case related to fQMCGGA. In addition, the gradi¬ 
ent correction (IDETGGA) to the IDETLDA form is needed 
to improve the performance for positron affinity calculations. 
Meanwhile, the IDETGGA correlation scheme makes dis¬ 
tinct improvement upon positron affinity calculations com¬ 
pared with the PHCGGA scheme which is similar to the cases 
of positron lifetime calculations of bulk materials. Neverthe¬ 
less, the best agreement between the calculated and experi¬ 
mental positron affinities is still given by the fQMCGGA & 
PZ-LDA approaches. 


V. CONCLUSION 

In this work, we probe the positron lifetimes and affini¬ 
ties utilizing two new positron-electron correlation schemes 
(IDETLDA & IDETGGA) based on three common electronic- 
structure calculation methods (ATSUP & ELAPW &PAW). 
Eirstly, we apply all approximation methods for electronic- 
structure and positron-state calculations to the cases of Si and 
Al, and give detailed analyses on the effects of these dif¬ 
ferent approaches. Especially, the difference between cal¬ 
culated lifetimes by using the self-consistent (ELAPW) and 
non-self-consistent (ATSUP) methods is clearly investigated 
in the view of positron and electron transfers. The well imple¬ 
mented PAW method with reconstruction of all-electron and 
full-potential, is found being able to give near-perfect agree¬ 
ment with the ELAPW method, which proves our calculations 
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TABLE V: Theoretical and experimental positron affinities (in unit of eV) based on four positron-electron correlation schemes and several 
electron structure calculation methods. The RMSDs between the theoretical and experimental positron affinities are also presented. Here, the 
PZ-LDA approach is labeled by PZ, and the PBE-LDA approach is labeled by PBE for short. 


IDFTGGA IDFTLDA PHCGGA fQMCGGA 



FLAPW 

PAW 

FLAPW 

PAW 

FLAPW 

PAW 

FLAPW 

PAW 

Exp. 


PBE 

PBE 

PZ 

PBE 

PBE 

PZ 

PBE 

PBE 

PZ 

PBE 

PBE 

PZ 


Si 

-6.481 

-6.478 

-6.683 

-6.884 

-6.881 

-7.070 

-6.728 

-6.726 

-6.926 

-6.182 

-6.179 

-6.373 

-6.2 

Al 

-4.497 

-4.504 

-4.683 

-4.624 

-4.631 

-4.813 

-4.641 

-4.648 

-4.828 

-3.981 

-3.988 

-4.169 

-4.1 

Fe 

-3.914 

-3.877 

-4.290 

-4.323 

-4.289 

-4.707 

-4.120 

-4.084 

-4.498 

-3.544 

-3.508 

-3.925 

-3.3 

Cu 

-4.381 

-4.437 

-4.932 

-4.875 

-4.933 

-5.435 

-4.614 

-4.671 

-5.168 

-4.073 

-4.130 

-4.630 

-4.3 

Nb 

-3.847 

-3.841 

-4.085 

-4.112 

-4.107 

-4.355 

-4.020 

-4.014 

-4.260 

-3.399 

-3.394 

-3.641 

-3.8 

Ag 

-5.147 

-5.083 

-5.577 

-5.670 

-5.615 

-6.109 

-5.398 

-5.337 

-5.831 

-4.875 

-4.817 

-5.310 

-5.2 

W 

-1.956 

-1.982 

-2.304 

-2.225 

-2.254 

-2.580 

-2.121 

-2.149 

-2.472 

-1.491 

-1.520 

-1.844 

-1.9 

Pb 

-5.954 

-5.936 

-6.305 

-6.328 

-6.305 

-6.683 

-6.186 

-6.166 

-6.538 

-5.622 

-5.601 

-5.977 

-6.1 

RMSD 

0.285 

0.283 

0.546 

0.570 

0.566 

0.899 

0.431 

0.427 

0.740 

0.314 

0.314 

0.272 

- 


are quite credible. While for ATSUP method, its competitive¬ 
ness against the FLAPW method is reduced within calcula¬ 
tions utilizing the best positron-electron correlation schemes 
(fQMCGGA). Then, we assess the two new positron-electron 
correlation schemes: the IDFTLDA form and the IDFTGGA 
form by using a reliable experimental data on the positron life¬ 
times and affinities of bulk materials. The gradient correction 
(IDFTGGA) to the IDFTLDA form introduced in this work is 
found necessary to promote the positron affinity and/or life¬ 
time calculations. Moreover, the IDFTGGA performs bet¬ 
ter than the PHCGGA scheme in both positron affinity and 
lifetime calculations. However, the best agreement between 
the calculated and experimental positron lifetimes/affinities is 
obtained by using the fQMCGGA positron-electron correla¬ 


tion scheme. Nevertheless, the new introduced gradient cor¬ 
rected correlation form (IDFTGGA) is currently competitive 
for positron lifetime and affinity calculations. 
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